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1. Introduction 

The present note is short version of our paper 1 where we have attempted a 
systematic comparison of the theoretical predictions of the /c t -factorization 
approach 2 ' 3 ' 4 with experimental data regarding the charm production pro- 
cesses at HERA. 

The production of open-flavored cc— pairs in ep-collisions is described 
in terms of the photon-gluon fusion mechanism. A generalization of the 
usual parton model to the fc t -factorization approach implies two essential 
steps. These are the introduction of unintegratcd gluon distributions and 
the modification of the gluon spin density matrix in the parton-level matrix 
elements. The hard scattering cross section for a boson gluon fusion process 
is written as a convolution of the partonic cross section &(x g , kt; 7*5* — ► qq) 
with the k t dependent (unintegratcd) gluon density A(x,kf 7 [i 2 ). 

The multidimensional integrations can be performed by means of 
Monte-Carlo technique either by using VEGAS 5 for the pure parton level 
calculations, or by using the full Monte Carlo event generator 
CASCADE 6 > 7 < 8 . 

Cross section calculations require an explicit representation of the fc t 
dependent (unintegratcd) gluon density A(x, /cf, /1 2 ). We have used three 
different representations, one (JB) coming from a leading-order perturba- 
tive solution of the BFKL equations 9 , the second set (JS) derived from a 
numerical solution of the CCFM equation 6 ' 7 and the third (KMR) from 
solution of a combination of the BFKL and DGLAP equations 10 . 

2. Numerical results and discussion 

A comparison between model predictions and data in principle has to be 
made on hadron level and only if it turns out that hadronization effects are 
small will a comparison to parton level predictions make sense. However, 
a full simulation even of the partonic final state, including the initial and 
final state QCD cascade needs a full Monte Carlo event generator. Such 
a Monte Carlo generator based on fc t -factorization and using explicitly off- 
shell matrix elements for the hard scattering process convoluted with k t - 
unintegrated gluon densities is presently only offered by the CASCADE 6 ' 7 ' 8 
program which uses the CCFM unintegrated gluon distribution. 

In the ref. 1 we systematically compared the predictions from the fct- 
factorization approach to published data on charm production at HERA. 
For this we use D* photo-production data from ZEUS 11 and D* produc- 
tion in deep inelastic scattering from both ZEUS 12 and HI 13 . First we 
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Figure 1. Relative differential cross sections l/oda/dx° bs (a) and l/odcr / d\cos 9*\ (b) 
with Q 2 < 1 GeV 2 , p v' > 3 GeV, \n D * < 1.5, 130 < W < 280 GeV, \r r ietl > 2 \ < 2.4, 
j^jet\,2 ^ ^ Mjj > 18 GeV and \fj\ < 0.7. The histograms are results with various MC 
simulations. 

calculate observables using a pure parton level calculation based on the 
matrix element calculation of BZ 14 including the Peterson fragmentation 
function 15 for the transition from the charm quark to the observed D* me- 
son. Then we compare the result with a full hadron level simulation using 
the Monte Carlo generator CASCADE with the matrix element of CCH 4 . 
We choose the JS unintegrated gluon for this comparison. 

Next we investigate on parton level different unintegrated gluon densi- 
ties. We study the sensitivity of the model predictions to the details of the 
unintegrated gluon density, the charm mass and the scale. We observed 1 
that the pt distribution of D* mesons both in photo-production and deep 
inelastic scattering is in general well described, both with the full hadron 
level simulation as implemented in CASCADE and also with the parton 
level calculation supplemented with the Peterson fragmentation function. 
We can thus conclude, that the pt distribution is only slightly dependent 
on the details of the charm fragmentation. 

We also consider the rapidity distribution of the produced D*. In photo- 
production and in DIS the differential cross section da/dr), where 77 is the 
pseudo-rapidity of the D* meson, is sensitive to the choice of the uninte- 
grated gluon distribution, We observed, that the parton level prediction 
including the Peterson fragmentation function is not able to describe the 
measurement over the full range of rj. The effect of a full hadron level sim- 
ulation is clearly visible as CASCADE provides a much better description 
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of the experimental data. 

Then we investigate the x 7 distribution, which is sensitive to the de- 
tails of the initial state cascade. We compare the predictions from a pure 
parton level calculation and a full event simulation of CASCADE with the 
measurements (Fig. 1(a)). We can conclude that the k t -factorization ap- 
proach effectively simulates heavy quark excitation and indeed the hardest 
p t emission comes frequently from a gluon in the initial state cascade 16 . 

Other interesting quantities are the dijet angular distributions of re- 
solved photon like events (x° BS < 0.75) compared to the direct-photon 
like events (x° BS > 0.75). In the fc t -factorization approach the angular 
distribution will be determined from the off-shell matrix element, which 
covers both scattering processes. Comparisions of the CASCADE results 
with the ZEUS experimental data for these angular distributions were done 
by S. Padhi 17 (Fig. 1(b)). 

In summary we have shown, that the k t - factorization approach can be 
consistently used to describe measurements of charm production at HERA, 
which are known to be not well reproduced in the collinear approach. We 
have also shown, that in fct-factorization, resolved photon like processes are 
effectively simulated including the proper angular distributions. 
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